Scott correction for large atoms and molecules in a 
self-generated magnetic field 



Laszlo Erdos * 
Institute of Mathematics, University of Munich 
Theresienstr. 39, D-80333 Munich, Germany 
S0ren Fournais ' 
Department of Mathematical Sciences, Aarhus University 
Ny Munkegade 118, DK-8000 Aarhus, Denmark 

and 

Jan Philip Solovej * 
Department of Mathematics, University of Copenhagen 
Universitetsparken 5, DK-2100 Copenhagen, Denmark 

May 2, 2011 

Abstract 

We consider a large neutral molecule with total nuclear charge Z in non-relativistic 
quantum mechanics with a self- generated classical electromagnetic field. To ensure sta- 
bility, we assume that Za 2 < kq for a sufficiently small kq, where a denotes the fine 
structure constant. We show that, in the simultaneous limit Z — > oo, a — > such that 
K = Za 2 is fixed, the ground state energy of the system is given by a two term expansion 
+ ci{k) Z 2 + o(Z 2 ) . The leading term is given by the non-magnetic Thomas-Fermi 
theory. Our result shows that the magnetic field affects only the second (so-called Scott) 
term in the expansion. 
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1 Introduction and the main result 

We introduce the molecular many-body Hamiltonian of N dynamical electrons and M static 
nuclei in three space dimensions. The electron coordinates are x\, x 2 , . . . , x^ £ M 3 , the nuclei 
are located at R = R 2 , ■ ■ ■ , Rm) £ ffi 3M . Let Z = (Zi, Z 2 , . . . , Z M ) denote the nuclear 
charges, Zj > 0, with total nuclear charge Z = Y2k=i Z^. We assume that the system is 
neutral, i.e. N = Z, in particular Z is integer. The electrons are subject to a self-generated 
magnetic field, B = V x A, where A £ iJ 1 (IR 3 ,IR 3 ) is a magnetic vector potential. The 
magnetic field energy is 

B 2 = ^ ( IV x A 



87ra 2 J R 3 8na 2 

where a is the fine structure constant. 

In the non-relativistic approximation, the kinetic energy operator of the j-th particle is 
given by the magnetic Schrodinger or the Pauli operator, 

TW(A) = (-iV Xj + A{xj)f, or T®(A) = [a- ■ (-iV Xj + A{x 3 ))]\ (1.1) 

depending on whether the particles are considered spinless or have spin--. Here cr = (<7i, a 2 , cr 3 ) 
is the vector of Pauli matrices. The Schrodinger operator acts on the space L 2 (IR 3 ), the Pauli 
operator acts on L 2 (IR 3 ,C 2 ). We will work with the Pauli operator, the treatment of the 
magnetic Schrodinger operator is simpler and we will only comment on the modifications 
needed. 

The electrostatic potential of the electrons is the difference of the nuclear attraction 

M 

k=l 

and electron-electron repulsion 

l<i<j<N lXt X ^ 

The total energy of the electrons is given by the Hamiltonian 

N 1 
H N (Z,R,A):=J2[T U) (A)-V(Z,R,x j )}+ J-^Tr (L2) 

j=l l<i<j<N 1 Jl 

This operator acts on the space of antisymmetric functions /\ x H, where % = L 2 (IR 3 ) ® C 2 is 
the single particle Hilbert space. 



\x — Ri 
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For a given vector potential A, the ground state energy of the electrons is given by 



E(Z,R,A) :=mfSpecH N (Z,R,A) (1.3) 

with N = Z = J2k Zk- By gauge invariance, E(Z, R, A) depends only on the magnetic field 
B = V x A. Considering the magnetic field dynamical, we focus on the absolute ground state 
energy of the system, that includes the field energy, 

£(Z,R,A) + -L f |v x A\ 2 , (1.4) 
Sua 1 J R 3 

and we will minimize over all vector potentials. 

Since we are interested in gauge invariant quantities (like energy, ground state density), 
we can always choose a divergence free gauge, V • A = 0. In this case, the field energy is given 
by 

|2 



|Vx^ = — — / |V<M ■ (1-5) 
87ra z J R3 87rcr J R3 

Since the magnetic energy will always be finite, we can also assume that A e L 6 (IR 3 ) (see 
Appendix of [FLL] for the existence of such a gauge), and we thus have 

V-A = 0, (f A 6 ) 1/3 <c[ |V®A| 2 = / |VxA| 2 (1.6) 



by the Sobolev inequality. 

We will call a vector potential A admissible if A e L 6 (M 3 ), V®Ae L 2 (M 3 ), and V • A = 0. 
For admissible vector potentials, the total energy is 

S(Z, R, A, a) := E(Z, R, A) + -i- f |V <8> A\ 2 , 

8na 2 J R3 

and the absolute ground state energy of the system is given by 

£ abs (Z,R,a) :=mf{f(Z,R,A,a)}, (1.7) 

where the infimum is taken over all admissible vector potentials A. 

Our units are h 2 (2me 2 )~ 1 for the length, 2me A h~ 2 for the energy and 2mech~ 1 for the 
magnetic vector potential, where m is the electron mass, e is the electron charge and H is 
the Planck constant. In these units, the only physical parameter that appears in the total 
Hamiltonian f ll.4p is the dimensionless fine structure constant a = e 2 (hc)~ l ~ We will 
assume that max^ Z^a 2 < k with some sufficiently small universal constant kq < 1 and we 
will investigate the simultaneous limit Z — > oo, a — > 0. 

The main result of this paper is: 
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Theorem 1.1. Fix M G N. Let z = (zi,z 2 , ■ ■ ■ , Zm) w ^h z x , z 2) • • • , -^m > 0, ^fc = 

and r = (r 1; r 2 , . . . , tm) G M 3Af u>z£/i min^ — > r min /or some r min > be given. With 
a positive real parameter Z > 0, define Z = (Z\, Z%, . . . , Zm), Z^ := Zzk, and R = Z~ 1/,3 r to 
be the charges and the locations of the nuclei. Then there exists a constant E TF (z,r) and a 
universal (independent of z, r and M), monotonically decreasing function S : (0, k,q] — > M. with 
some universal k > and with lim K ^ + = | such that as Z = X^feLi %k — * oo, a — > 
with maxfc 87rZfca 2 < k q , we have 



M 



£ abs (Z, R, a) = Z 7/3 E TF {z, r) + 2Z 2 ^ ^(St^o: 2 ) + o{Z 2 ). (1.8) 



fe=i 



Remark 1.2. The constant E TF (z,r) is given by the Thomas-Fermi theory for nonrelativistic 
molecules without magnetic field, see below. The factor 2 in the Scott term is due to the spin 
degeneracy. 

Remark 1.3. Our current proof does not provide an effective error term in (11. 8p although we 
conjecture that o(Z 2 ) can be replaced with 0(Z 2 ~ V ) with some rj > where the constant in the 
error term depends only on r m i n and M. 

Remark 1.4. From the monotonicity of £ , a b s (Z, R, a) in a, it is obvious that the function S 
is monotonically decreasing. It is known that there is a finite critical constant n cr > such 
that the system is unstable if max^ Z^a 2 > n cr , in particular the restriction max/; Z^oi 2 < 
in the theorem is necessary but our threshold k is not optimal. We conjecture that S(k) is 
a strictly decreasing function, but unfortunately our proof does not provide this statement. In 
fact, we even cannot exclude the possibility that S(k) is constant (=1/8) for all k up to the 
critical value n cr beyond which it is minus infinity. 

Remark 1.5. An energy expansion similar to (II. 8p was derived in a relativistic model (without 
magnetic fields) in ISSSf . In the relativistic case one must also consider simultaneous limits 
Z — > oo, a 0. In this case, however, it is the combination Z^a which must remain bounded 
in contrast to Theorem \l.l\ where a bound on Z^a 2 is required. In the atomic case, M = 1, 
an alternative proof of the relativistic energy asymptotics was given in IFSWltf . 

The scaling can be understood from Thomas- Fermi theory which we recall briefly \LS\ [L] . 
Let < g(x) G L 5 / 3 (IR 3 ) n L 1 (M 3 ) then the Thomas- Fermi energy functional is defined as 

£ TF ( ) : = ^(3vr 2 ) 2 / 3 / g 5/3 - [ V(z, r, x)g(x)dx + D(g) 
5 Jr3 Jm.3 



with 



1 Jr3 J r s \x - y\ 
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(the coefficient in front of the kinetic energy takes into account the spin degeneracy). It is 
well known that the variational problem 

. M 

mi{s TF ( 6 ) : QEL^R^nL^iR 3 ), /e = X>*} (1.9) 

J k=i 

has a unique, strictly positive minimizer, called the Thomas-Fermi density and denoted by 
g TF (x) = g TF (r,z,x). The value of the minimum, E TF (z,r) := £ TF (£> TF ), is called the 
Thomas-Fermi energy. The function 

V TF (z,r,x) := V(z,r,x) - g TF * \x\-~ 1 (1.10) 

is called the Thomas-Fermi potential] it is strictly positive and it solves the Thomas-Fermi 
equation 

K7 = (3- 2 ) 2/3 K] 2/3 - (1-11) 
Sometimes we will use the notation Vj F (x) instead of V TF (z,r,x) and likewise for £>J F . The 
key quantities in the Thomas-Fermi theory have the following scaling behavior: 

V TF (z, r, x) = /T 4 V TF (/i 3 z, h'h, h~ l x) (1.12) 
£ TF (z, r, x) = h~ 6 Q TF (ti 3 z, frh, h~ l x) 
E TF (z,r) = /i- 7 E TF (/i 3 z,/i" 1 r) 

for any h > 0. We also note that the Thomas- Fermi energy defined as the minimal value of 
(jl.9p can also be given by the phase-space integral of the classical symbol with the Thomas- 
Fermi potential, i.e. 

E TF (z, r) := 2-1- / [p 2 - V^( q )]_dpd q - D(£ F ), (1.13) 

where [a]_ = — min{a, 0} denotes the negative part of a real number or a selfadjoint operator. 
The factor 2 in fll . 13j) accounts for the spin degeneracy. 

Suppose we replace the many-body electrostatic potential in Hn(Z, R, A) by its mean field 
approximation as 

N 1 

3=1 l<i<j<N 1 n 

Then the absolute ground state energy £ , a b s (Z, R, a) will be approximated by 

inf {Tr [T(A) - _ + -L £ | V ® A\ 2 - D(eif R )} (1.14) 
« Z 4 / 3 inf {Tr [T h (A) - V™]_ + -^hT* |V ® A| 2 } - Z 7 / 3 D(,J F ), 
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by using the scaling property (11.121) with the choice h = Z l ^ 3 and introducing the notation 



T h {A) = [<r-{-ihV + A)} 2 . 

The infimum in (11.141) is taken over all admissible vector potentials A, but in fact taking 
infimum over slightly different sets yields the same result, see a more detailed discussion in 
Appendix A of [EFST] . 

Note that the trace in (11.141) is computed in the spinor space L 2 (IR 3 , C 2 ). Any operator M. 
acting on L 2 (M 3 ) will be naturally identified with acting on the spinor space L 2 (M 3 ) £*DC 2 

and we will not distinguish between M. and M. <8> I in the notation. This will not cause any 
confusion but we need to keep in mind that Tr l 2 (r 3 )®c 2 -M = 2Tr l 2 (r 3 )-M- Unless indicated 
otherwise, the traces in this paper are computed on L 2 (M 3 ) ® C 2 . 

We thus reduced the problem to a semiclassical analysis of the Pauli operator with a self- 
generated magnetic field. The leading term in the asymptotic expansion in negative powers 
of h = Z~ x l 3 for the one particle energy in HI . 14[) is given by the Weyl term 

^{Tr^-^^/^lV^I 2 } 

= h ~"wv I ^ ~ V™(q)]-dpdq + o{h~ 3 ) (1.15) 

as long as k := 8irZa 2 < k,q with some sufficiently small fixed positive k,q (in fact, it is 
sufficient that max^ 8irZka 2 < k,q). Notice that the leading term does not depend on k. Via 
( II. 13ft this produces the leading term asymptotics Z 7 ^ 3 E TF (z, r) + o(Z 7 ^ 3 ) for the many body 
ground state energy. 

A fundamental result of Lieb and Simon |LS] (see also [L]) rigorously justified this heuristics 
without magnetic field and with an effective error term. Thus they proved the leading term 
asymptotics in ( II. 8p for the ground state energy of large atoms and molecules without magnetic 
field. The next order term, known as the Scott correction, is of order Z 2 . For the non-magnetic 
case it is explicitly given by 




(1.16) 



fc=i 



(the additional factor 2 is due to the spin degeneracy) and it was rigorously proved for atoms in 
[H ] ISWlj and for molecules in [IS] (see also |SW2[ I5W3| ISSj). The next term in the expansion 
of order Z 5 ^ 3 was obtained in |FSj . 

It was established in |ES3j that the inclusion of a self-generated magnetic field does not 
change the leading term asymptotics in (II. 8p . The main theorem in this paper, Theorem 11.11 
shows that the effect of the magnetic field appears in the second (Scott) term in the asymptotic 
expansion. 

The proof has two main steps. First we reduce the interacting many-body problem to an 
effective one-body semiclassical problem as it was described in (j!.14p . This step will be done 



6 



rigorously in Section [3J The second step is an accurate second order semiclassical asymptotics 
for the magnetic problem. More precisely, we will show the following more accurate version 
of dUTSJ): 

inf{Tr [T h (A)-V?J]_ + ^ /jV®A| 2 } (1.17) 

9 r - 

= h ~ 3 7^-s / & - C(?)]-dpdg + 2/^ 2 £ z^(8vrZ fc a 2 ) + (/i~ 2 ), 

with k = 8irZa 2 and under the assumption that max^ 87rZfcQ; 2 < n with some k that is 
sufficiently small. Together with (OS]) . (OIj) and /i = Z" 1 / 3 , (071) will yield (Ojl . 

The precise second order semiclassical expansion for 

infjTr [T h (A)-V] + / |V®A| 2 ) (1.18) 

A I J Ktl 2 J R3 ) 

is of interest itself for a general potential V and with k < Kq. Under general conditions on 
V the leading term is given by the classical Weyl term as in (|1.15|) . A local version of this 
statement was proven in Theorem 1.3 of |ES3] (this theorem contains only the lower bound, 
the upper bound is trivial). The global version was given in Theorem 2.2 of |EFSlj . where 
the condition k < k could even be relaxed to k = o(/i~ 1 ). 

The subleading term in the expansion for (11. 18ft depends on the singularity structure of 
the potential. For V £ we proved in Theorem 1.1 of |EFS2] that the Weyl asymptotics 
holds with an error 0(h~ 2+ri ). The main ingredient was a local semiclassical asymptotics that 
we recall in Theorem 12.21 

Using this result for smooth potentials and a multiscale analysis, we will show in Section [2] 
that for potentials with Coulomb singularities, a non-zero second order term arises from the 
non-semiclassical effects of the innermost shells of the Hydrogen-like atoms. The precise form 
of this term will not be as explicit as in the non- magnetic case, (jl.lfip . but it will be given by 
a universal function S which we will describe in Theorem 12.41 and prove in Section HI 



2 Semiclassical results up to second order 

In this section we are interested in noninteracting fermions, each is subject to the one-body 
Hamiltonian H = Th(A) — V. The total ground state energy of the system is given by 

5(A) :=Tr [T h (A)-V)_ + ^- [ | V x A\ 2 = Tr [T h (A) - V] _ + / |V®A| 2 , (2.1) 

Kh' J R 3 Kh* J R 3 

where k > is a parameter, and where the last equality uses that V • A = which can be 
assumed by gauge invariance. 
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We will assume that the potential has a multiscale structure. Intuitively, this means that 
there exist two scaling functions, /, £ : M 3 — > M + such that for any u G M 3 , within the ball 
B u (£(u)) centered at u with radius £(u), the size of V is of order f 2 (u) and V varies on scale 
£{u). Moreover, we also require that the continuous family of balls B u (£{u)) supports a regular 
partition of unity. The following lemma states this condition precisely. This statement was 
proved in Theorem 22 of [SS] with an explicit construction^. 

We will use the notation B x {r) for the ball of radius r and center at x and if x = 0, we 
use B{r) = B (r). 

Lemma 2.1 ([SS] Theorem 22]). Fix a cutoff function ip G C^°(M 3 ) supported in the unit ball 
B(l) satisfying J> 2 = 1. Let £ : R 3 ->■ (0, 1] be a C 1 function with ||W||oo < 1- Let J(x,u) 
be the Jacobian of the map u (-> (x — u)/£(u) and we define 



ip u (x) = if) 



x — u 

Ku) 



^/J(x,u)£{u) 3/2 . 



Then, for all x G 



ip u (x) 2 £(u)- 3 du = 1, 



\n\ 



+ n 3 , 



and for all multi-indices n G N 3 we have 

Halloo < cj{u)-w, 

where C n depends on the derivatives of ip but is independent of u. 
We will require that the potential satisfies 

\d n V{u)\ <C n f(u) 2 £(u)-^ 



(2.2) 

(2.3) 
□ 

(2.4) 



for all n G N 3 uniformly in u in some domain Q C M. 3 . In applications, Q will exclude an 
/z-neighborhood of the core of the Coulomb potentials. For brevity, we will often use £ u = £{u) 
and /„ = f(u). 

Inserting the partition of unity (12.21) . by IMS formula and reallocation of the localization 
error, we have 



£{A) =Tr 
>Tr 

> 



— U u [T h (A)-V]ip u -K 2 \Vilj u \ 



du 



ij u [T h (A) -V - Ch 2 \Vip u \ 2 }ij u 



+ 
+ 



Kh 2 
1 

^2 



|V® A\ 2 
IV® A\ 2 



^£(a,v u + m 



(2.5) 



1 Multiscaling was introduced in semiclassical problems in [IS] (see also |Sob[ ISobl] ) 
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where 

V+:=V + Ch 2 \V^ u \ 2 

and 

£(A,V+M:=Tr\i/> u [T h (A)-V+Wu\ + 4? / ^ ® A ? ■ 

L J — KH J R3 

In (12.51) we used Tr [j O u du]^ > J Tr [0 M ]_du for any continuous family of operators O u . 
We will assume that V£ satisfies the same bound (12.41) as V, i.e. 

h < CfJ u . (2.6) 

For the Coulomb-like singularity, V(x) = l/\x\, we will choose the £(u) = c\u\ with some 

c < 1 and f u = £u 1/2 so that (J23D holds. With this choice (J3ID holds for \ u \ > Ch<2 ^ so the 
multiscale analysis will work apart from a very small neighborhood of the nucleus. 

Next we recall our local semiclassical result from |EFS2] on a model problem living in the 
ball B(£) of radius I > 0. 

Theorem 2.2 (Semiclassical asymptotics [EFS21 Theorem 3.1]). There exist universal con- 
stants no G N and e > such that the following is satisfied. Let K ,f,£,ho > and let 
k < K f- 2 £-\ Let ip e C™(R 3 ) with supp^ e B(£) and let V G C°°(B{£)) be a real valued 
potential satisfying 

\d n i/j\ < C n £~ H , \d n V\ < C n f 2 £~ H (2.7) 
for every multiindex n with \n\ < Uq. Then 



inf (Tr [i)H h {A)^}_ + -\- I \ V ® A\ 2 ) - 2(2vr/i)- 3 / ^{qf [p 2 - V(q)] _dqdp 

< Ch~ 2 + £ f- e £ 2 - £ (2. 



for any h < h f£. The constant C depends only on k q , h and on the constants C n , in (12. 7p . 
The factor 2 in front of the semiclassical term accounts for the spin and it is present only for 
the Pauli case. □ 

Remark. By variation of k, we obtain from Theorem 12.21 the following estimate 

B{21) 



Kh 2 

for (near) minimizing vector potentials A 



A\ 2 < Ch~ 2+£ p- £ £ 2 - £ , (2.9) 



The following result from |EFS2] can be viewed as a partial converse to (12. 9 p as it estimates 
the semiclassical error in terms of the magnetic field. Note that the assumption in ( I2.10p below 
is much weaker than (12. 9p . 
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Theorem 2.3 ( |EFS2} Theorem 3.2]). Let the assumptions be as in Theorem \2.2\ and assume 
that A satisfies the bound 



>B(2£) 

Then, with e from Theorem \2.2\ we have 



[ |V® A\ 2 <Ch~ 2 f£ 3 . 
JB(2e) 



(2.10) 



Ch~ 2 f£ 3/2 { [ |V® A| 

Jb(2£) 



, V2 

+ Ch~ l ft 



> Tr [tpH h (A)iP}_ - 2(27r/i)- 3 / iP(q) 2 [p 2 - V(q)] dqdp 

> Ch- 2+£ f- £ £ 2 ~ £ - Ch~ 2 f 2 £ [ |V ® A| 2 , 

where the constants may depend on ho and kq and on the constant in (12.101) . 



(2.11) 



Very close to the nuclei (at a distance h 3 ~ Z~ l ) the semiclassical approximation is no 
longer valid due to the local Coulomb singularity and the energy is given by a specific function 
S that describes the contribution of the innermost electrons. Since the nuclei are separated 
on the much larger semiclassical scale h, this effect is additive for different nuclei. Thus S 
can be defined via the Hydrogen atom. The following theorem defines S and gives some of its 
properties. It will be proven in Section HI 



Theorem 2.4. Let <p e C£°(R 3 ) be a cutoff function with supp0 C B(l), = 1 on 5(1/2) 
and such that <\> := (1 — (f) 2 ) 1 / 2 G C°°(IR 3 ). Define 4>r(x) = 4>(x/R) for some R > 0. There is a 
universal critical constant n cr > such that for any k < K cr and for any fixed < < (2k;) -1 
the following limit exists, it is finite and it depends only on k ( and it is independent of the 
choice of (ft and (3 satisfying the conditions listed above): 



lim 



inf \ Tr 

A 



R[ T h=1 (A)- — 



- [ |V® A\ 2 + f3 [ \V®A\ 2 

K JB{R/i) J B(2R)\B{R/A) 



2(2tt)- 



4>Uq) 



? ■ 



dpdq 



-:2S{k). (2.12) 



The limit is denoted by 2S(n) following the convention in the literature to indicate explicitly 
the factor 2 that accounts for the spin degeneracy. The function S(k) is defined on [0, n cr ), 
it is decreasing and 5"(0) = 1/8, corresponding to the coefficient of the non-magnetic Scott 
correction. 
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Furthermore, there exists a sequence {Ar} C i? 1 (IR 3 ) with supp A R C B(R/A) and such 
that 



lim 



Tr 



T h=1 (A R ) - i.) 



1 

K J B(R/A) 



\V ®A R \ 



- 2(2tt) 



P _ rr 



dpdg 



25(«). (2.13) 



This way of introducing a Scott correction when one cannot calculate its explicit value 
was introduced in |Sob] and was used later in |SS| ISSS] . In the case of the relativistic Scott 
correction an alternative method of characterizing the Scott term was given in [FSWll IFSW2] . 

Equivalently, one could define S(k) via another limiting procedure, which may look more 
canonical: 

Lemma 2.5. For the function S(k) defined in ( 12.12ft we have 



lim 



inf i Tr 

A 



T h=1 (A) - -1- + v 
\x\ 



[ \V®A\ 2 } 



2(2tt) 



1 

\q\ 



dpdq 



2S{k). (2.14) 



In the main arguments of this paper we will always use the definition ( 12. 121) and only 
in Section EJ with the help of the main semiclassical theorem, Theorem 12.71 we will prove 
Lemma 12.51 

Now we formulate the precise second order semiclassical asymptotics for a potential with 
Coulomb like singularities and with certain scaling properties that are satisfied by the Thomas- 
Fermi potential ^ TF . We first specify the properties of V TF that are used in the proof. 



For any given r = (r l5 r 2 , . . . , r M ) G M 3 and z = (zi,z 2 ,.-., z M ) e 



set 



mm \ru — rA 



and 



d(x) := mm{\x — rk\ : k = 1, 2, . . . , M} 



f(x) :=mm{d(x)- 1/2 ,d(x)- 2 }. 



(2.15) 



If M = 1, then we set r min = oo. We say that a potential V is of Thomas-Fermi type if it 
satisfies the following two properties: 
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(i) There exists \x > such that for any multiindex a with \a\ < n (where the universal 
constant n Q is given in Theorem 12. 2p we have 

d:[V(z,r,x)+fi}\ < C a f(x) 2 d(x)^ a \ (2.16) 

where the constants C a depend only a, M and max^ z^, 

(ii) For \x- r k \ < r min /2, we have 

-C< V(z, r, x) - < Cr~l + C (2.17) 

\X T k\ 

where C depends on z. 
Then we have 

Theorem 2.6 ( |SSt Theorem 7]). The Thomas-Fermi potential V = V TF satisfies the condi- 
tions (EHI and fl2Tfjl . □ 

With these ingredients in Section [5] we will prove the following second order semiclassical 
asymptotics: 

Theorem 2.7 (Semiclassical asymptotics with Scott term). There exist universal constants 
n G N and k > such that the following is satisfied. Let V be a real valued potential 
satisfying (12.161) and ( 12.171) . Then 



lim 

h->0 



inf hi [T h (A)-V}^ + ^- [ |V®A 

A \ Kh z J R 3 

-2(2nh)~ 3 / [p 2 -V(q)]_dqdp-2h- 2 Vz^(z« 



(2.18) 



for any < k < kq . 

Moreover, there exist some e > and ho > and there exist an admissible vector potential 
A and a density matrix 7, whose density satisfies 

£? 7 < / [Vf! 2 + Ch~ 2+£ (2.19) 



3tt 2 



and 



d(q^ - (?>ix 2 )- l h- z [Vfl 2 ^ < Ch~ 5+£ (2.20) 
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Tr [T h (A) - V\<y + / IV ® ^4| 2 (2.21] 



/or any < h < h , such that 

Kh 2 „ 

<2(2nh)- 3 f J [p 2 -V(q)]_dqdp + 2h- 2 J2 z l S ( z kK) + o(h- 2 ). 

k=l 

The constants C in the right hand side of these estimates depend only on k , h and on the 
constants in (12.161) and (12.171) . The factor 2 in front of the semiclassical term accounts for 
the spin and it is present only for the Pauli case. 

Convention: All integrals, unless specified otherwise, are on M 3 . 



3 Proof of the Main Theorem 11.1 



In this section we complete the proof of our main theorem. 
3.1 Lower bound 

The first step is to reduce the many body problem to a one body problem. We will use the 
following Lemma whose proof relies on the Lieb- Oxford inequality |LO] . 

Lemma 3.1. \ES3\ Lemma 4-3] There is a universal constant Cq > such that for any 
G Af Cn R3 ) ® C 2 with ||tf|| 2 = I, for any non-negative function p : R 3 — > R with 
D(p, p) < oo, for any compactly supported and admissible A and for any 5 > we have 

N _ i , , r 

|2 



\5J2t^(a)+j2) — - — i]y) + c °f \^ xA 

1 i=1 iKj \ x i x j\ j 1 JW> 



N 



> -D(p,p) + h,J2(P* kil" 1 )*) -CS^N. □ 



i=i 



Thus, applying this lemma for g = g^, the density of we get 
{V,H N (Z,IL,A)V) 



> (^,y"Ul-5)T^(A)-V(Z,R,x J )]^) + D(g li ,)-C5^ 1 N -C [ |V x A\ 2 

N „ 

>h,Y / [(l-S)T^\A)-V TF (Z,R,x j )]^)-D(gl F R )-CS- 1 N-C / |V x A\ 

> Tr [(1 - 6)T(A) - V?l] - D(gl F R ) - CS^N - cj |V x A\ 2 , 

JM. 3 
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where in the second step we used the definition fll.lOp and that D(q^ — q^b) > 0. We now 
add the field energy (11.51) and absorb the —Co J | V x A\ 2 term at the expense of factor (1 — 5) 
by using 

8na 2 87TQ! 2 

as long as 5 ^> Z^ 1 and Za 2 is bounded. Now we use the scaling properties (11.121) of the 
Thomas-Fermi theory with h := Z^ 1 ^ 3 and k := SnZa 2 to get 



£ abs (Z,R,a 

>z 7/3 



h 3 (l - 25) inf (Tr [T h (A) - V^]_ + -L J | V ® A| 2 ) - £>(£ 
^^inf (Tr [T fc (A) - 2l/f r F ]_ + J \V ® A 



C^iV 



>z 7/3 



M 



P 2 - V™(q) dqdp + 2hJ2 4S(z k K) - o(h) - D{ B , 



k=i 



> z 1 ' 3 



(2tt) 3 
C5Z 7/3 - CS^N 

M 

E™ + 2hJ2 z 2 k S(8nZ k a 2 ) - C5 - o(h) 



z,r I 



k=l 



- C5^N 



M 



>E^ + 2Z 2 J24S(SnZ 



kOt 



\Z 2 ) - CZ 2 ^ 6 . 



(3.1) 



k=l 



From the second to the third line we used ( 12. 18ft from Theorem 12.71 twice. In the main term 
the 1 — 25 prefactor can be easily removed for a lower bound since the term inf^ ( • • • ) is non- 
positive. In the error term we used that inf^ (...) is bounded by 0(h~ 3 ) since f[V^J] 5 / 2 < C 
with a constant depending only on M (see (12. 16ft ). Then we used (11. 13ft and finally the scaling 
relation fll . 12j) for the energy. In the last step we also inserted the optimal 5 = Z~ 5 ^ 6 and 
used N = Z. This completes the proof of the lower bound in Theorem ll.il 



3.2 Upper bound 

By Lieb's variational principle |L2] and neglecting the exchange term, after a rescaling by 
h = Z~ x l 3 the energy of the particles (II. 3p can be estimated by 



E(Z, R, A) < Z^ 3 (Tr [T h (A) - V(z, r, -)]7 + ZD(Z~ 1 g 1 ) 

= z^ (tt [T h (A) - v™] 1 + zd(z~% - ell) - ZD (C)) ( 3 - 2 ) 
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for any density matrix 7 on L 2 (IR 3 ,C 2 ) with density g y (x) = Tr C 27(x,x) and with Tr 7 = 
J g(x) < N = Z. Adding the field energy, we have 

E abs (Z,R,a)<Z 4 / 3 (Tr [T h {A) - VjJ^ + -^— 2 h' 2 [ |V ® A\ 2 

+ ZD(Z-i ej -e™)-ZD(e%j) (3.3) 



for any admissible vector potential A. 

Let k = 8irZa 2 , set V = V z T r F and choose an admissible vector potential A and a density 
matrix 7 according to (I2.19p . (I2.20p and (I2.2ip . In particular 



Tr 7 



= < z(\ + cz- l/3 ~ £/3 ) 



81 



using that 



3tt 2 



TFn3/2 



by (fTTTj) and the constraint / g = J2 k z k = 1 in CT> - We thus defme 7 = (l + CZ" 1 ^^ 3 )" 1 ^ 
so that the constraint Tr 7 < Z is satisfied. Moreover, from (12.201) we have 

D(Z-ig^-g™)<CZ-^-^. 
From the triangle inequality for \/D we obtain 

D{Z- 1 g 1 - gll) < C(l + CZ-^l 3 fD(Z- x g~, - g™) + CZ^^D^™) < CZ^^/ 3 
using that -D^JJ) < C and s < 1. In summary, we obtained 

£ abs (Z,R,«) <^ 4/3 (Tr [UA)-VjJ] 1 + ^- 2 h^ J \V®A\ 2 -ZD(gZ)+CZ^^) 



for the vector potential ^4 from Theorem 12.71 Since (12.211) holds for 7 and the energy of 7 and 
7 differ by a factor (1 + CZ -1 / 3- ^ 3 ), we obtain, after the usual rescaling, 



£ abs (Z,R,a) 

<2 7/3 



2( 27 r)- 3 || b 2 -C F (^)]-d^-^Ka 



2Z 2 J> 2 S(^) + o(Z 2 ). 



k=l 



Using the identity f ll . 13j) . we thus obtain the upper bound in (II. 8p which completes the proof 
of Theorem O □ 
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4 The Scott term 

In this section we give the proof of Theorem 

Proof of Theorem \2.4\ First we prove the existence of the limit (12. 1 2 p . Define, for R > (and 
with k, (3 as in the statement of the theorem) 



>R 



IV® A\ 



IV ® A\ 



K JB(R/4) 
2(27t)" 3 



(4.1) 



B(2R)\B{R/A) 



(?) 



V 



and 



1 



S(R,k,P) :=-ini£ R , K> ?(A). 



(4.2) 



We will prove later in (14. lip that this infimum is not minus infinity. 
Step 1: A-priori upper bound. 

Upon inserting A = we get S(R, k, (3) < S(R, 0, oo). Since we know that limR-j.oo S(R, 0, oo) = 
S(0) exists |SSS| Lemma 4.3 with a = 0], we obtain that S(R,k,/3) is bounded from above 
uniformly in R. I.e., there exists a constant K (independent of R, k and (3) such that 



S(R,k,P) < S(R,0,oo) < K . 



(4.3) 



Step 2: Lower bound and semiclassics. 

Consider r < R/8. Let <p r satisfy = 1 an d define W r = |V0 r | 2 + |V0 r | 2 and <fi r> R = 4>R((>r- 

Split 

r + <f>r,R(T h=1 (A) - ^)4> r ,R 



>> R {T h=l (A) - > <Pr{T h=l (A) ~ r 

\x\ \x 



■ (T h= i(A) - |— | 



r ,fl(T h=1 (A) - — 



- W r 



>r,R 



■ (4-4) 



The second term will be estimated by borrowing a small part of the field energy. We will use 
that the local semiclassical result Theorem 12.21 holds with any positive constant in front of 
the field energy. However, this regime has to be treated with multiscaling. The proof of the 
following lemma is postponed to the end of this section. 

Lemma 4.1. For any r > and 5 > and for any r, R satisfying r < r < R/8 we have 



inf < Tr 

A 



r , R (T h=1 (A) 



1 

\x\ 



>r,R 



5 

I B(2R)\B(r/4) 

> 2(2tt)~ 3 



|V® A\ 





r 2 1 1 




P ~ VI 


i \x\ J 



(4.5) 



dxdp — Cs t1 . r £ ^ 2 , 
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where e > is the exponent obtained from Theorem \2.2\ and the constant in the error term 
depends only on r and S. 

Therefore, combining (14.41) and (14. 5p . we have 



Tr 



^ R (T h=l (A)-— )<j> R 



dxdp 

|V ® A\ 2 - C s , ro r- e/2 

< B(2R)\B(r/4) 

Hence, combining the semi classical integrals, for any ro < r < R/8 we get 



> Tr 
-5 



(T h=l (A)- — -W r 
\x\ 





\ 2 1 1 


l <t> 2 r,R( X ) 


P ~ Tl 


L \x\ J 



K JB(R/4) 



|V® A\ 2 + (3 



(4.6) 



IV® A\ 



B(2R)\B(R/A) 



-5 



|V® A\' 



B(2R)\B(r/A) 



2(2vr)~ 3 / 



2 1 1 

\x\ 



dxdp - C«5 iro r 6/2 . 



(4.7) 



In the following, we will fix an ro > 0. 
Step 3: A-priori lower bound. 

We fix r = ro- On the ball of radius ro we use that if n < k* where k* is a small universal 
constant, then 



Tr 



and 



<j) ro (T h=1 (A) - 



-2(2vr) 



1 



1 

2k 



\V®A\ 2 > -K x 



B(2r ) 





r 2 1 1 






L i J 



dxdp > —K%, 



(4.* 



(4.9) 



where the constants only depend on r . The estimate (14. 9 p follows by simple integration. The 
other estimate (14.81) is a consequence of (the proof of) [ES3| Lemma 2.1] in the case Z — 1. 
Inserting ffl~8l) and ffl~9]) in (jUTJ), we get 

Sr,kAA)>(^--S) [ |V® A| 2 + (/3-<5) / |V® A| 2 -^(r ,5) (4.10) 

^ K JB(R/2) J B(2R)\B(R/2) 



for any k < K* and with a finite constant i^(ro, 5) depending only on ro and 5. Choosing 
5 = (3 < (2k)" 1 , this proves in particular that 



S(R, k, (3) > — oo, 



(4.11) 
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and with the choice 5 = (2k,*) 1 we get 

1 S( j R,k,(2k)- 1 ) > -C(k*) (4.12) 
for some finite function C(k*) < oo depending only on r$ and n* as long as K < k*. 
Step 4: Bound on the magnetic energy. 

From (14.101) we also obtain that if A is such that it yields a better energy on the large ball 
than no-magnetic field, i.e. 

S R ,^{A)<£ R ^{A^Q), 

and k < k*, then 

/ |V® A\ 2 <C{r ,/3) (4.13) 

JB(2R) 

with some constant depending only on tq and f3. This follows from (I4.10p (taking 5 = <5o(/3) := 
i min((2K*) _1 , f3)) because then we have 

W) 



[ |V ® A| 2 < ^, Ki/3 (A = 0) + K(r , 5 (P)), 

Jb(2R) 



'B(2R) 

and £^ iKi( a(A = 0) is known to have a finite limit as R — > oo, independent of k, /3 \SSS\ 
Lemma 4.3]. 

Step 5: Bound on the localization error. 

We will now remove the localisation error W r from (14. 7\\ using the bound on the field energy. 
By the variational principle and W r < Cr~ 2 , we can estimate 



Tr 



■(T h=1 (A) 



- W r 



\.r\ 



> Tr 



CV~ 2 Tr l ( _ OO)0 )(^(T ft=1 (A) - - W r )(j) r y (4.14) 



We have, with g 
Tr 1 



-{|i|<r}; 



1 



(p r (T h=1 (A) - — - W r )(f) r j < Tr l(-oo,o)| 



{p + Af 

• 1 { -oo,0)[(P + A) 2 - 



1 

L \x\ 



— - \B\ - Cr 

\x\ 

\B\ + Cr~ 



9, 



Here we used the fact that we consider the strictly negative eigenvalues to get the last in- 
equality. By the CLR estimate, we therefore have 



Tr 1 



(-oo,0) 



1 



)<cf 


- 1 


- \B\ + Cr~ 2 


H\x\<r} 


-\x\ 





3/2 



dx 



{\x\<r} 



\X 



3/2 



CC"\r 3 / 2 + r 3 ^ 



\B\ 3 ' 2 + Cr- 3 

|2\ 3/4 



dx 



B 



{\x\<r} 



) 
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where we used the Holder inequality to get the last inequality. Using the uniform bound (I4.13P 
on the field energy, we can therefore control the last term in (14.141) for any r > r as 



r- 2 Tr l(-oo,o) ((j>r (T h=1 (A) - -L - W r )(j> r ) < C ro r^ 2 



(4.15) 



Thus 



Tr 



,{T h=1 {A)- — -W r )(t> r 



x\ 



> Tr 



^.(T h=1 (A)- — 
\x\ 



C r j- X l\ (4.16) 



Step 6: Monotonicity of the energy in the radius. 

Combining (14. 7p . (I4.16P and the definition (14. ip with R replaced with r and replaced 
with P' we get 



1 

+ 



-6/2 



+P / -P' / -s 

K J B(R/4)\B{r/A) J B(2R)\B(R/&) J B{2r)\B{r/4) J B{2R)\B{r /A) 



>g r ^,(A) + ^ [ 

z JE 



iv® A\ 2 -a 



B(2R)\B(r/A) 



So,ro' 



-e/2 



|V® A\ 2 
(4.17) 



for any (3', (3 < 1/(2k) and if we choose 5 = 5q{k, (3) := | min{(2fi;)~ 1 , (3} and recall that 
r < R/8. Notice that we even saved a part of the field energy in B(2R) \ B(r/4). Taking 
infimum over all A, we have 



S(R, K,P)>S{r, K,(3')-C So , ro r~ £/2 



(4.18) 



which shows that R — > S(R, k, /3) is essentially an increasing function with a uniform upper 
bound (14. 3p . hence it has a limit. To be more precise, define S(k, (3) := limsup^^ S(R, k, (3), 
then for any rj > there is a sufficiently large r = r(rj) > tq such that Cs . ro r~ £ ^ 2 < rj/2 and 
S(r,K,P) > S(k,/3)-t]/2. Then ffl~TM implies that S(R,k,P) > S(k, /3)-t] for any R > Sr(r)). 
Together with the definition of S(n, (3) this means that S(n, /3) = lim^oo S(R, n, /3). 

Step 7: Independence of the limit of (3. 

Suppose that (3 < f3' < 1/(2k). Then clearly S(R, k, (3) < S(R, k, /?'). Furthermore, by taking 
first the limit R — > oo in (I4.18p . then the limit r — > oo we obtain that 

lim S(R,k,(3) > lim S(R,k,P'). 

R—^oo R—>oo 



So we get that the limit S(k, (3) is indeed independent of (3 < 1/(2k). Moreover, from (14.181) 
it follows that 

S(K)>S(r,K,P)-Cs , ro r- £/2 (4.19) 
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for any r > r , (3 < 1/(2k) and 5 (k, (3) : = |min{(2fi;) Combining this bound with 

(14. 121) . we obtain in particular that 

inf S(k) > C(K*,r ) > -oo (4.20) 

K<K* 

for some constant depending only on n* and r . Together with the upper bound H4.3|) this 
shows that S(k) is a bounded function for k e (0, k*]. 

The fact that £(0) = |, i.e. the non-magnetic case, has been proven before, see e.g. 
Lemma 4.3 in [SSSj with the choice a = 0. (Note that in |SSSj S(0) = \ is stated but there 
the kinetic energy was — |A, while our non-magnetic kinetic energy is —A which accounts for 
the apparent discrepancy.) 

It remains to prove that one can obtain S(n) by considering only vector potentials with 
small support. 

Step 8: Improved bound on field energy. 

Let A R be an (almost) minimizer of the variational problem (14.21) . i.e. 

2S(R,k,(3)>£ r ^(A r )-R^ 2 . 

Using fT4TT7D with the choice r = R/8, (3' = (3 and estimating £ R / BtKt p(A R ) > 2S(R/8, k, (3), 
we get 

S(R, k, (3) > S(R/8, k, /3) + § / |V <g> A R \ 2 - Cs ,r R~ e/2 (4.21) 

Z J B(2R)\B(R/32) 

for any R > 8r . Now letting R — > oo, we conclude that 

lim - / |V® A R \ 2 = 0. (4.22) 



R-toc 2 



B (2 R)\B(R/ '32) 



Step 9: Upper bound on S(R,k,/3). 

Fix r = R/8 and k, (3. By the definition of S(r, k, (3), there exists a vector potential A r such 
that 

2S(r,K,P) >E r ^p{A r )-r-\ (4.23) 

and we can assume that 

A r = (4.24) 



Ibi 



' B(2r)\B(r/4) 

by adding a constant to A r if necessary. Finally, by (I4.22p we have 



e(r) := / |V® A r \ 2 = o(l), (4.25) 

' B(2r)\B{r/32) 



as t — V oo. 
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Furthermore, there exists a density matrix 7 r such that 
2S{r,K,p) >Tt (p r -f r (f) r (T h=1 (A r ) - — ) +- / \\7®A r \ 2 

\ \X\J K JB(r/4) 



(4.26) 



P 



IV® A 



'B(r/4) 
|2 o/o_\-3 



2(2tt) 



We define A' 



' B(2r)\B(r/4) 

^rA, then A' r = A r on the support of <p T . Moreover 



2 1 1 

V i r 
9 



dpdg — 2r . 



B(2r)\B(r/4) 



|V<8>4.| 2 < Cr" 2 

< c 



V ® A r 



B(2r)\B(r/A) 



(4.27) 



f B(2r)\B(r/4) 

|V<8> A r \ 2 < Ce{r) 

' B(2r)\B(r/A) 

with a universal constant C. Here we used the Poincare inequality on the ring B(2r) \B(r/4) 
which holds with a universal constant since the width of the ring is comparable with its radius. 
Thus, we have from (I4.26[) 



2S(r,K,(3)>Ti ( j )rlr( f ) JT h=1 (A' r )-^-)+- [ \V ® A' r \ 2 

V \X\/ K Jfl( r /4) 



(4.28) 



2(2vr) 



■(q) 



? ■ 



dpdq — C{r 1 + e(r)). 



We use this A' r as a trial vector potential for S(R, k, (3), i.e. we have 



2S(R,k,P) <Tr (j> R [T h=1 {A! r 



\x\ 



- 2(2tt) 



(9) 



P 2 - 



+ 



K JB(2r) 

dpdq. 



/ |2 



(4.29) 



The field energy integral (with coefficient k~ x in (14. ip ) can be restricted to B(2r) and the 
second field energy integral is absent since A' r is supported on B(2r). 

Now we construct a suitable trial density matrix 7. It will have the form 



7 := (firlAr + (j>r,Rj& 



r,R, 



with < 7 < 1 to be chosen below. Inserting 7 into ( I4.29P and using the support properties 
of the cutoff functions, we obtain 

2S(R, K,P) <£ R , K A A 'r) 

<Tr (f> rlr( f> r (T h=1 (A' r ) - A) + Tr ^rI^rIt^A^ - y 



V® A'X -2{2n) 



B{2r) 



9 



dpdg 



<2S(r, k, /3) + C(e(r) + r" 1 ) + A(t), 



(4.30) 
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where 



A (7) := Tr <p r ,RH ( ( ) r,R\Th=i(A' r 



\x\ 



2(2vr) 



p 



1 1 

I Q 



dpdq. (4.3i; 



We used (14231) . flOH|l and (14271) in the last step of ( IPOD . 

We will use Theorem 12.31 and multiscaling to prove that we may choose 7 such that 

A(t) ^^v/iW + r" 1 / 2 ). (4.32) 

Inserting (fl~32D in ( H~30j) and recalling that e(r) ^ by (1425]) we obtain (1243]) (with the 
choice Ar = A' r ). Here we used the choice that r = R/8 and the previous result that 
S(k) = lim^oo S(R, k, (3) exists and is independent of j3. 

It thus remains to prove (I4.32p for a suitable choice of 7. Let ip u be given as in Lemma [241 
with the choice £(u) = £ u := i5q\Ao + u2 - Define 



7 := 



where 
and 



V :={x : r/3 < \x\ < R} 

1 



7« 



-(-00,0] 



lpu4>r,R{T h=1 (A' r 



\x\ 



Clearly < 7 < 1 using (TO) . 

Inserting this choice of 7 into (14.291) and using ( 12.21) . we obtain 

A(t) = / -^{Tr [^.^(TW^) - o)^^ 



- 2(2 



7T 



i-3 



V 



1 

l<?l J 



dpdg 



(4.33) 



Notice that we could restrict the du integration to V since otherwise ip u 4>r,R vanishes due to 
the support properties of these functions. 
We will use Theorem 12.31 for each u with 

h = i, f~ 2 = e u . 



We have 



IV® A'\ 2 < C 1 



B{u,2l u ) 



V® A'X < Ce(r) 



B(2R)\B(r/4) 



by ( I4.27p . Here we used that suppA^. C B{2r) to get the last inequality. So (12. 101) is satisfied 
and we get from ( 12. lip and (14.331) that 



A(7)<C^|{e(r) 1 / 2 + C 1/2 }<^(e(r) 1 
for r > r . This finishes the proof of (14 . 3 2 p . 



/2 _l r -l/2 



(4.34) 
□ 
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Proof of Lemma \4.1\ We choose £{u) 
functions and define the ring 



h\frl + u 2 and f u 



-1/2 



100 



for the scaling 



V :={x : r/3 < \x\ < R} 



which supports <p r ,R- Inserting the partition of unity ( 12. 2 p and reallocating the localization 
error we get 



Tr 



(j)r, R (Tl{A) - — - W r )(j) r , R 



Tr 



du 



v 



2 ±2 



.R 



du 



i*u<l>r,R\Tl{A) 



X 



C(W r + \V^ u \ 2 ))<P r ,R^ 



Notice that we could restrict the du integration to V since otherwise ip u 4> r ,R vanishes due to the 
support properties of these functions. We also used that Tr [J O u dfi(u)}_ > J Tr [O u ]_d/i(-u) 
for any continuous family of operators O u and for any measure [i. We can also reallocate the 
field energy as 



S B\ 



> B(2R)\B(r/4) 

with some positive universal constant c. Thus 
1 

\x 



■ I ^ f |V® A\ 2 



Tr 



I |V<8> A\ 2 > 

B{2R)\B(r/A) 



du 



r,R.{Tl(A) - — - W r )i 

where we define 

E r , R (A,U^ u ) := Tr i) u <p r)R (Ti(A) - U)(p rtR ^ u 

for any potential U, and in the last step we used that 

|1 + c(w r + < + j) =■ v + {x). 



S r , R (A,V + ^ u ), (4.35) 



: / |V 

B u {2i u ) 



A\ 



This inequality holds for any u G V by the support and scaling properties of W r and ip u and 
by the estimate W r < Cr~ 2 l(r/2 < \x\ < r) < Cr^ 1 /\x\. 

It is easy to see that ip = ip u 4>r,R and V = V + satisfy the condition (12. 7p in Theorem [2^2] on 
the ball B u (2£ u ) (the theorem was formulated for balls about the origin but it clearly holds for 
balls with different center). We can choose h > 1 and r sufficiently large so that c5 > Rq 1 
and 1 < h f u £ u = ^/io(l + u 2 ) 1 ^ are satisfied for all u G V and r > r . Thus, Theorem 12.21 
with h = 1 gives 

£r, R (A,V + M > 2(2ny 3 J J [(^ r , R )(x)} 2 [p 2 - V + (x)}„dxdp-C 5 , ro £Z £ / 2 , 
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where the constant C^ ro is independent of u but it depends on 5 and on r (via the choice of 
k and h ). 

Inserting this bound into (I4.35p . we get 

r,Il(Tl(A) 



Tr 



W r 



\.~)'\ 



r,R 



6 



V <g) A\ 



B(2R)\B(r/A) 



>2(2vr)" 3 1 1 ^ R {x) 2 [p 2 -V + {x)}Mp-C & , ro [ 



where in the first term we extended the du integration from V to R 3 and we used (12. 2\\ . The 
last integral is bounded by C^r -2 / 2 , uniformly in R. 

Finally, we can remove the localization errors from V + , since 



<j)r, R ( x ) 2 ([p 2 - v + (x)}^ - [ P 2 - i^rvWdp 



f 1 r/ C\5/2 

/, , ( 1 + 7) - 1 

J x \>r 2 \ x \ ' 1 v ' x 



< c 

'|x|>r/2 \ X \ 

< Cr" 1 ' 2 



dx 



(4.36) 



with a constant independent of R. This error can be absorbed into the other error as e < 1. 
Thus we get (14. 5 p and have proved Lemma 14.11 □ 



5 Semiclassics with Scott term 

In this section we prove Theorem 12.71 The guiding principle follows the similar proofs in 
[SS[ ISSSj . We first divide the space into three regions. The first region consists of disjoint 
balls of radius r ~ h about the nuclei. Here we will use the Scott asymptotics as described 
in Theorem 12.41 The second region is far away from the nuclei, at a distance R > h~ x . The 
contribution of this regime will be estimated by a simple Lieb-Thirring inequality. Finally, in 
the third intermediate region, we will use an argument similar to Lemma [4.11 which relied on 
the multiscale decomposition and Theorem 12.21 on each domain. The precise decomposition 
is the following. 

Choose two localization functions 9± £ C 1 (1R) with the properties that < 9± < 1, 
B 2 _ + 9\ = 1, moreover 6L(t) = 1 for t < 1/2 and 6L(t) = for t > 1. Recall that d(x) denotes 
the distance of x to the nearest nucleus and r min is the minimal distance among the nuclei. 

For any r < r min /A and R > r min we set 

± (*):=e ± (^), * ±W :=* ± (^). 

Note that 

M 

(f)-(x) = y^ y 9 r>k (x), with 6 r>k (x) = 9-(\x - r k \/r). 

k=l 
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Assuming h is small enough, we will choose 



Ch~ 1 ' 2 if [L = 
CR^ if /i ^ 



(5-1] 



with some small £ > 0. Here the //-dependent constant i? M > is chosen such that — V(x) > 
for d(x) > RJ2. Clearly + ®\ = 1, 0?. + = 1 and 

This latter partition of unity corresponds to the three regions we described above. The 
localization errors will be affordable in all regimes. Near the nuclei, the localization errors 
of order h 2 r~ 2 = 0(h 2 ^) is more than one order of magnitude smaller than the size of the 
potential in this regime, which is \x — r^" 1 ~ h -1 , and the localization error relative to the 
potential becomes even weaker further away from the nuclei. The localization error far away 
is of order h 2 R~ 2 = h 3 will be negligible both in L 5//2 and L A norms as necessary for the 
magnetic Lieb-Thirring inequality that we recall for convenience: 

Theorem 5.1. ILLSf There exist a universal constant C such that for the semiclassical Pauli 
operator T h (A) - V with a potential V £ L 5/2 (1R 3 ) n L 4 (M 3 ) and magnetic field B = V x A £ 
L 2 dR 3 ) we have 



Tv[T h (A)-V]_<Ch~ 3 f [V]*l 2 + C(h- 2 J Wf\J [V][ 



1/4 



(5.2) 



5.1 Lower bound 

For any e £ (0, 1/4), we define 

%{A) : = Tr [e r ,k(T h (A) - V - Ch 2 r - 2 )6. 
1 



r.fc 



+ 1 



e\ 1 



2/ nh 2 



|V® A\' 



Br k (r/4) 



+ 



AkK 2 



IV® A\ 



(5.3) 



'B rfc (2r)\B rfc (r/4) 

Using the IMS localization and the fact that the balls B rk (2r) are disjoint since r < r m j n /4, 
we have 



Tr [T h (A)-V}- + -^ [ |V®A| : 

Kll J R3 



M 
k=l 



+ Tr \<5> + {T h {A) - V - Ch 2 W R )$- 



+ 



2nh 2 



IV® A\ 



+ Tr 



$-<f>+(T h (A) -V- Ch 2 W r , R )<t>+®. 



(5.4) 
IV® A\ 2 



K 3 \U fc ^ fc (r/4) 
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with some positive universal constant C and with Wr := |V$_| 2 + [V^+j 2 and W r ^ 
|V0_| 2 + |V0+| 2 + |V$„| 2 + |V$+| 2 . 

First line in (15. 4p 

Fix k and recall from (I2.17P that 

-V{x) > 



\x - r k 



~ Cr min — C 



on the support of 6 r ^- Thus we can write 

r k (A)>T k {1 \A) + T h {2) (A) + T^(A), 

where 



T k (1) (A) :=Tr 6 r , k ((l - 2e)T h (A) - (1 - 2e> 

L \x-r k \ 

l-2e 



+ 



1 -2e 



S r , (2r) 



8k/i 2 



|V® A| 



B r . (2r)\B_ (r/4) 



T fc (2) (A) :=Tr 6 r , k (eT h {A) - Ch 2 r~ 2 - Cr^ n - C)d r , k 



Anh 2 



IV® A\ 



B r J2r) 



T k (3 \A) :=Tr 6 rtk (eT h (A) - 2e-^-)6 r 



x - r k 



Anh 2 



IV® A\ 2 . 



B r , (2r) 



After pulling out the common (1 — 2e) factor, shifting r k to the origin and rescaling, we obtain 
from (1212)1 . with /3 = (8k)" 1 , that 



infT fc (1) (A)>(l-2e) 



2(2tt/i)" 



fl 2 

"Vfc 



V 



\Q — rk\ 



dqdp + 2h z k S(z k K 



o{h~ 



Here we made use of the fact that after rescaling the variable R in (12.121) becomes rh^ 1 = h~^, 
so the limit h — > is equivalent to the limit i? — > oo in (12.121) . 
The term in the square bracket can be bounded by 



2(2tt/i)" 



2 

r.k 



V 



z k 



\Q - r k\ i 



dqdp + 2h 2 z 2 S{z k n) 
dq 



\g-r k \<r \q-r k \ 5/2 



+ Ch~ 



<Ch~ 3 

<Ch- 3 r^ 2 + Ch~ 2 , 
using z k < 1 and S is bounded (if k < k,q is sufficiently small). Thus 



(5.5) 



inf T k {1) (A) > 



2(2irh)' 



9 



P 



Zk 



\q - r k \i 



dqdp + 2/i 2 z k S{z k K 
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as long as er 1//2 <C h. 

The term t£ 2 \A) is estimated by the magnetic Lieb-Thirring inequality ( 15. 2 p which we 
need in the following form: 

Lemma 5.2. Let <p G C^°(]R 3 ) be a cutoff function with supp <p C B{1), = 1 on 5(1/2). 
Define 4>i(x) = <j)(x/£) for some I > and /e£ f2 := supp<p£. Then for any vector potential A 
we have 

Tr [MT h (A) - V)<f> e ]_ + [ |V ® A| 2 > -C7T 3 / Kf /2 - C/T 3 / (5.6) 



with some universal constant C. 



Proof. Let <p e (x) := <fi(x/£), where e C~(M 3 ) is supported in -8(3/2) and <p = 1 on -B(l). 
Set (A) := |_B(2£)| _1 J B / 2 ^) ^ anc ^ ^ := ~ (A))4>£- By a gauge transformation and by the 
fact that A = A — (A) on the support of <pg, we have 

Tr [MTh(A) - V)<f>i]_ = Tr [fa(T h (A) - V)<f>z]_. 

By the magnetic Lieb-Thirring inequality (I5.2p we have 

Tr [MTk(A) - V)cf> e ]_ > -Ch~ 3 Jjl' 2 - c(h~ 2 J |V x I| 2 ) 3/4 ( jvf) (5.7) 
By supp A C B(2£) and by the Poincare inequality 



\\7xA\ z = \W®A\ 2 <C \V®A\ 2 + Cr 2 \A-{A)\ 2 <C \W®A\ 2 . 

JR 3 J B{2£) Jb(2£) J B{2£) 

Combining this with ( 15. 7ft we obtain ( 15. 6p . □ 

We return to the estimate of 7^ {A), and shifting to the origin and using ( 15.61) we 
obtain 

7f (A) > - Ceh~ 3 I (V 1 [h 2 r~ 2 + + l]f 2 -Cej (V 1 [/> 2 r~ 2 + + l] ) 4 
> - Ce- z/2 h- z [h^r- 2 + r 3 ) - Ce~ 3 (h s r- 5 + r 3 ) > -Ce~ 3/2 h-^ - Ce~ 3 h 3 -^ 

with a constant C depending on k and r min . This error term is of order o(h 2 ) as long as 
e > h and £ < 1/6, which we will assume from now on. 

(3) 

Finally, for the term 7^. , after shifting, rescaling and using Zk < 1 

fr 2/i _1 i 1 r 

Tr 9 d (T h=1 (A) - —)9 d +— |V O A| 

Fl J q:h ' n J B(2d) 
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with d := rh^ 1 = and 9d{x) = 9-(x/d). Now we will use the "running energy scale" 
argument as in the proof of Lemma 2.1 of [ES3j . where 2/i -1 plays the role of Z in Lemma 
2.1 of |ES3j . the localization errors in (2.14) of |ES3] are not present, and the key condition in 
|ES3] that Za 2 is sufficiently small translates into k being sufficiently small. Using the final 
result (3.23) from |ES3] , and noting that the d~ 2 term there was due to the localization error 
that is not present now, we obtain 



inf7; (3) (A) 



> _ C£h -5/2 d X/2 



(5. 



which if o(h~ 2 ) provided e <C h^ +1 ^ 2 . These constraints, together with the previous conditions 
er 1 / 2 <C h, e > h and £ < 1/6 leave plenty of room, choosing for example £ = ^ and e = h 3 / 4 . 
In summary, for the first line in (15.41) we proved that 



A I 



i\l 



fc=i 



k=l 



2(2tt/z)" 



ti 2 

a r,k 



V 



Zk 

\q - r k\ 



dqdp + 2h z k S(zkH 



+ o(h~ 



(5.9) 



Second line in ( 15.41) 



By the magnetic Lieb-Thirring inequality from ( 15. 6ft we have 

e 



Tr 



$ + (T h {A)-V -Ch 2 W R )$. 



2nh 2 



|V<g> A[ 



> - Ch~ 3 [ [V + Ch 2 W R fl 2 - Ck 3 e~ 3 [ [V + Ch 
Jn+ Jo. + 



2tj/ 14 

+ ' 



(5.10) 



where fl + := supp$ + . The contribution of the W R terms is negligible using ||Wr||oo < CR 2 : 
and that its support has a volume CR 3 : 



h~ 3 / [h 2 W R f 2 + k 3 s~ 3 / [h 2 W R ] 4 < Ch 2 R~ 2 + CK 3 e~ 3 h 8 R- 5 < 



Ch 3 if/i = 
Ch 2 if ^ 

with a constant that may depend on \i. The positive part of the potential [V(x)]+ is zero for 
/i 7^ in Q + and it can be estimated by f(x) 2 < d(x)~ 4 according to (" 12 . 1 6[) if \i = 0, so 



h-*\y + tffi + tfe-*\y + & ) < 



C(h~ 3 R~ 7 + e~ 3 R~ 13 ) if ^ = 
if 7^ 



with a constant depending on Ko and M. With the choice of ( 15. ip and recalling e > h, the 
lower bound on the second line in ( 15. 4p thus vanishes as h — > 0. 
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Third line in 



This estimate will be very similar to the proof of Lemma 14. 1[ so we will skip some details 
here. We choose £(u) = l u := \Jr 2 + d(u) 2 and f u = minj^ 1 ^ 2 ,^ 2 } for the scaling 
functions and define the regime 



Q:={x : |z| < 2R, \x - r k \ > r/3, k = 1, 2, . . . , M} 



(5.11) 



which supports $-0+. Inserting the partition of unity (12. 2p and reallocating the localization 
error, we have 



Tr 



$_0 + (T,(A) - V - Ch 2 W r , R )<P + <$>- 
du 



8nh 2 



IV® A| 



K 3 \U fe S rfc (r/4) 



> 



Q 



where we defined 

£(A,U,jP u ) :=Tr 



Th{A) — U)(f) + &-ip u 



(5.12) 



+ 



Kill 2 



IV® A\ 



B u {2£u) 



for any potential U and with some sufficiently small positive universal constant c > 0, and we 
defined 

V+(x) := V(x) + Ch\W r . R {x) + |V^(x)| 2 ) 

on the support of ip u . We recall from (12.161) that V(x) < Cf(x) 2 < Cf 2 on the support of ip u , 
since f(x), defined in (I2.15p . is comparable with f u . Thus 



\V+(x)\ < Cf u + Ch% 2 < Cfl < Cf(xf 



on supp i/j u , 



where we distinguished the case \i = and /i ^ 0. In the latter case \u\ is bounded (depending 
on fi) and thus f u is bounded from below. We also used that h < C£ u f u = C mm{il/ 2 , i" 1 } , 
which holds since i u is between a constant multiple of r and R if u G Q. Similar estimate 
holds for the derivatives of V^ + , i.e. the main condition (12.71) of Theorem 12.21 is satisfied for 
with scaling functions I = £ u and / = f u . The other condition, k < Kof^i" 1 , is trivially 
satisfied. Applying Theorem 12.21 we get from (15.121) that 



inf 

A 



Tr 



> 



$„0 + (T h ,(A) - V - Ch 2 W T ^+<$>_ 



- + 8nh 2 J R 


|V® A\ 




3 \U k Br k (r/4) 



^h(2nh)- 3 j J [(^$_0 + )(g)] 2 ^-K + (9)]„dgdp-C'/ i - 2+ Vr e e 
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The second term is 0(h~ 2+e ) since the integral is finite even after extending to IR 3 from 
Q. In the first term we use that the localization errors in V£ are bounded by Ch 2 l~ 2 < 
Ch 2 [d(x) + r]~ 2 and are supported in a ball of radius CR, and thus 



h~ 3 / (^$_0 + ) 2 K + ]f <(1 +e)h~3 (^^mT + Ce-^h- 3 



I\x\<cr +rJ 
<h~ 3 J (^^ + ) 2 [V}f + Ceh~ 3 + Ce~ 3 / 2 h 2 R}' 2 



5/2 

dx 
(5.13) 



since \V\ < f 2 E L 5 ^ 2 . Choosing e = h 1+< > with a small ( > 0, and recalling that R < Ch 1 / 2 , 
we obtain that the two error terms are of order /i~ 24 ^, which even after the J Q £~ 3 du integration 
is o(h- 2 ). 

In the main term, we perform the du integration and use (\2.2\i to obtain 

Third line of Q > 2(27ih)- 3 J J [($_0+)(g)] 2 [p 2 - V(q)] _dqdp - o{h~ 2 ). 



Collecting the estimates of all three terms in H5.4[) and using the properties of the cutoff 
functions, we have 



Tr [T h (A) — V}--\ — ^— [ \V®A\ 2 

Ktl J R 3 

> 2(2nhy 3 I! [p 2 - V(q)] dqdp + 2h~ 2 4S(z k n) 
+ 2(27ih)- 3 J j <$>+{q) 2 [p 2 - V(q)]_dqdp 



M 



M 



2 Z k 
p - 



\Q ~ r k\ 



k=l 

The middle term in the r.h.s. in absolute value is bounded by 
Ch~ 3 



p 2 -V{q) \dqdp- o(h~ 2 ). (5.14) 



x\>R/2 



,5/2 < I ch~ 3 j lxl>R/2 d( x y w < ch- 3 Ri < ch'/ 2 if fx = o 

^ if fl ^ 0. 



The last term in (15.141) . also in absolute value, is bounded by 

M 



k=l 



~k 



\Q - r k \ 



5/2 



[V(q)} 



5/2 



dq < Ch~ 3 I |g|~ 3/2 dg = Ch~ 3 r 3/2 

W\<r 



by using (12.171) . where C depends on M and r m i n . With our choice of r = h 1 iq, this error 
term is also negligible. This completes the proof of the lower bound in Theorem 12.71 
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5.2 Upper bound 

We again set £ u = >/r 2 + d{u) 2 and consider the appropriate cutoff functions ifj u from f 12 . 2 j) . 
We construct a trial density matrix of the form 

7 = ^ &r,klk&r,k + / -^-0+^7^n0+, (5.15) 

fc=l J s *« 

where Q was defined in (15.1 ip and 7& and 7 M are density matrices to be determined below. 
Since 

we obtain from (12.21) that 7 is a density matrix. 



5.2.1 Trial density near the nuclei 

To construct 7^, we fix some 77 > and we notice that from the last part of Theorem 12.41 for 
any unsealed cutoff function 0, there exists some R{rf) such that for any R > R(rj) there is 
a vector potential Aq, supported in B(Rq/A), and there is a density matrix 7 such that 



Tr 



^4> Ro (T h=1 (A ) - 



+ 
2(2?r 



1 

Z k K 
3 



IV® A 



4,(9) 



? ■ 



We can also assume that 



Tr 



7<tao( r /i=i(4)) 



1 

l.rl 



dpdg < 2^(2 fc K) + 77. (5.16) 



(5.17) 



1 /2 

by noticing that the semiclassical integral is of order R which dominates over the 2S(z k K) +t] 
term for sufficiently large Rq. 

Fixing an appropriately large Rq, 7& and A k are now obtained from 7 and Aq by shifting 
and rescaling, such that 



Tr 



j k e r>k [T h (A k ) - 

- 2{2nhy 



Zk 



x - r k \ 
-3 



6 r j 



1 



kK 2 



P 



|V® A k \ 

z k 



\Q - r fe| J 



dpdq < 2h- 2 z 2 k S(z k K) + Ch~ 2 r] (5.18) 



with r = RQh 2 z k l . Here we used that the unsealed cutoff function can be chosen to be 
4>(x) = 8-(d(x)) so that after rescaling and shift became 9 r>k . We choose R = /i -1- ^ 
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such that r = h 1 * and clearly Rq > R(rj) is satisfied in the limit as h — > 0. We also 
remark that A k is supported in B r u(r k ) which are disjoint balls for different fc's. Defining 
A := Y^k=i Ak, we have A = A k in the support of 6 r ^. Thus summing up (I5.18p . and using 
that the replacement of z k \q — r k \~ l with V(q) in the semiclassical integral term is negligible 
(see the estimate of the last term in (15.141) ). we have 



M 



^Tr \ lk e r jT h (A) 



~k 



k=i 



\x - r k \ 



Kh 2 



IV® A\ 



M 



< 2(2tt/i)- 3 / <j) 2 _{q) p 2 -V(q) dpdq + 2/T 2 ^ z 2 k S{z k K) + Chr\ (5.19) 



fc=i 



T k (A) from (15.81) . we have for any L > 



Now we establish some properties of the density Qk(x) : = Jk( x , x )- Similarly to the estimate 



inf < Tr 

A 



L\Th{A) — 



\x-y\ 



nil 2 



IV® A\ 2 \ > -Ch~ 3 L 1/2 



uniformly for any y 6 M 3 if k is sufficiently small. In particular, for any density matrix 7 with 
density g 1 we have 

2 SUp / 02 (x) ^L dx < C </,-3 jL l/2 + Tr70LT/i(A)0L + _L / |V®A| 2 

for any vector potential A and thus 



sup 



\x - y\ 



dx < Ch~ 3 L 1/2 + inf <^ Tr 



l4>L (r h{A) i)J + _L| 



IV® A| 



Applying this bound to 7 constructed above (with L = R , h = 1 and A = A ) and using 
(j5.17|) . we have 



sup / ^{x)^-dx<CR { 



,1/2 



x - y\ 

which, after rescaling and shifting 7 to 7^ amounts to 



sup / 2 ik (x) 

yeK 3 J 



QiA x ) 
\x - y\ 



dx < Ch~ 3 r 1/2 . 



(5.20) 



In particular, 



/ el k {x) B , k {x)dx<Ch~^' 2 



(5.21) 
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by choosing y = r k and using that \x — r k \ < Cr on the support of 6 r ^. 

Combining ( I5.2ip with (I2.17P we see that z k \x — r k \~ l can be replaced with V in the l.h.s. 
of (I5.19P at an error of order h~ 3 r 3/2 = o(h~ 2 ) and thus we have 



^Tr [ lk e rjk (r h (A)-vy rjk 

< 2(2nh)- 3 



k=i 



Kh 2 



IV® A| 



M 



(q)p 2 -V(q) dpdq + 2h~ 2 J24 S ( z k^ + Ch- 2 7]. (5.22) 



k=i 



Moreover, it follows from (I5.20p and (15.21 j) that the density 9 2 ^ k g lk of the density matrix 
9 r ,klk6r,k satisfies 



D (9l,kQik 



e r Ax) 2 e r , k {y) 26lk{x)6lk{v) A ~ 



x - y\ 



-dxdy < Ch~ 6 r 2 . 



(5.23) 



5.2.2 Trial density away from the nuclei 



Now we construct 7„ for any u G Q. Since supp A k C -B r / 4 (rfc) and (j)+(x) is supported at 
d(x) > r/2, thus on the support of <p + we have A = 0. Therefore it is sufficient to construct 
a non-magnetic trial state 7„ within each ball B u (£ u ) which supports ip u . This was achieved 
in Corollary 15 of |SS] and we just quote the relevant upper bound (we formulate it particles 
with spin, this accounts for an additional factor 2 compared with |SSj ): 

Proposition 5.3. JSS\ Corollary 15] Let \ £ C*o(^ 3 ) be supported in Bi with some i > and 
let V G C 3 (Bi) be a real potential. Assume that for any multiindex n G N 3 with \n\ < 7 we 
have 

\\d n X \\oo < CJ-W, HcTVIU < C n f 2 r^ (5.24) 
with some constant f . Then there exists a density matrix 7 such that 

Tr [ 1X (-h 2 A-V) X ] <2(2irhy 3 J J x(<?)V ~ V(q)}^dpdq + Ch~ 3+6/5 f 3+4/5 £ 3 ~ 6/5 , 

where C depends only on the constants in (I5.24p . Moreover, the density q 1 {x) of ^ satisfies 
Ql {x) - 2(27rh)- 3 io 3 [V(x)}f < C/l -3 + 9/io / 3-9/io r 9/iO ; (5>25) 

for almost all x G Bi, and 



X 2 g y -2(2nh)- 3 u 3 / X 2 [V] 



i3/2 



< Ch~ 3+6/5 , 



(5.26) 



where u 3 = 4n/3 is the volume of the unit ball. 
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We will apply this Proposition for each ball B u (£ u ) and with \ := <p + ip u . From (I2.3P and 
(I2.16P it is easy to see that the conditions ( I5.24p are satisfied with our choice of £ = i u = 
^Q-y/r 2 + d(u) 2 and /„ = f{u) = min{ci(u)~ 1 / 2 , d(u)~ 2 } as given in (I2.15p . The density matrix 
thus constructed in Proposition 15.31 will be denoted by 7^ and this will be the density matrix 
in (I5.15p . With this choice we have 



Tr 



f 

'-h 2 A-V) J —<j) + ^ u j u ij; u (f) + 

< 2(27r/ i )- 3 J ^ J J MrfMrflp 2 - v(q)].d P d q + ch-^ 5 J ^/„ 3+4/5 e 6/5 

< 2(2nh)- 3 J J MQ) 2 [p 2 -V(q)Udpd q + Ch^ | s ^/„ 3+4/5 C 6/5 - (5-27) 
The error term can easily be computed by using that £ u ~ d(u) for u G Q as 

f ^ f u +A 'Hl-^ < c [ dtuy^-^du + c f d{uy^i^du 

JQ "u Jr/3<d(u)<l Jl<d(u)<2R 

< Cr- 1/W + C. (5.28) 

Since r > h, we have / i 1 /s r - 1 / 1 o < h 1 ^ 10 , so the error term in ( 15.27ft is o{h~ 2 ). 

We now use 7 from (15 . 1 5[) and A = A k constructed in Section 15.2.11 to complete the 
proof of the upper bound (I2.2ip in Theorem 12.71 Combining (I5.22p and (I5.27P and recalling 
that A and 0+ are supported disjointly, we have 

Tr[T h (A)-V} 1 + ^ J \V®A\ 2 

r r M 

< 2(2vr/i)- 3 J J [P 2 ~ V(q)] _dqdp + 2h~ 2 £ z 2 k S(z k K) + o(h~ 2 ) + Ch~\ 



where we also used that <p 2 + + = 1. Finally, letting h — > first and then rj — > 0, we obtain 
the upper bound in (l2~i8|) and (12721"]) . 

To complete the proof of Theorem 12.71 it remains to prove (I2.19P and (12.201) . From (15.151) 

we have ^ 

k=i 

Using ( 15 .211) and 0+ < 1 we have 
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The first term is of order h 3 / 2 ( 1+ £) = h 2+7 / 20 ; hence negligible. In the second term we use 
(Ell and to get 



du 

Q 



3vr 2 /l 3 In 11 



du 



Q 



du 



(5.29) 



Since both logarithms are of order | logh\, we obtained (I2.19p . 

To prove (I2.20p . we note that WD(gj satisfies the triangle inequality, thus we have 



du 



2/2 



Q~tu 



3tt 2 /i 3 



[V] 



3/2 



(5.30) 



+ CD 



1 



3tt 2 /i 3 



[V] 



3/2 



du 



The last term is smaller than 0(h 5+e ) with some small e by using (I5.23p . For the second 
term in (15.301) we note that by (12.21) 



du 

Q ^\ 



ij u {xf<P + {xf - 1 < |0+(x) 2 -l| + 



/ ^u{x) 2 <C(l(d(x)<r) + l(d(x)>R)) 



and thus 



D 



3tt 2 /i 3 



[V] 



3/2 



73-^+ - 1 



D[[V] 3 _l 2 l(d{x) <r))+D[\V]% z l(d(x) >R) 



i3/2. 



By the Hardy-Littlewood-Sobolev inequality we have D(g) < C\\q\\qi 5 for any real function g, 
therefore this error is bounded by 



CK 



[V]f + 



d(x)<r 



d(x)>R 



Using that V(x) is essentially Zk\x — rjt| 1 near the nuclei and V(a;) ~ \x\ 4 for large x (see 
( 12.16P and (I2.17P ). we easily obtain 



D 



1 



,3lT 2 h 3 

which is smaller than 0(h~ 5+£ 



[V] 



3/2 



^^-l])<^- 6 (r 2 + i?- 7 ) 
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Finally, we estimate the first term on the r.h.s. of (15.30 p . By the Hardy-Littlewood-Sobolev 
inequality and (I5.25P 



3ir 2 h 3 



[V] 



3/2 



1/2 



<c 



, du 



du 2 2 



3n 2 h 3 



[V] 



3/2 



6/5 



1 



8ju 



3n 2 h 3 



[V] 



3/2 



<ch^ J p 2 u 1/10 i 8 J 5 



6/5 



<ch~ 2 -^ 

as the last integral is bounded. In the second line we estimated 



(5.31; 



1 



3n 2 h 3 



[V] 



3/2 



6/5 



by using (I5.25P and (12. 3p . This completes the proof of (I2.20p and the proof of Theorem 12. 71 □ 



6 Equivalence of the two definitions of S(k) 

Proof of Lemma \2.5[ Setting h = z/ 1 / 2 , after a change of variables we have 



(27T) 



p - — + v 

m 



dpdq = v 



(2nhY 



P -T1 + 1 



dpdq. 



Similarly, by a simple rescaling, x — > v 1 x we get 



inf i Tr 

A 



Ti{A). 



-v 



X 



+- / |V®Ar \ = i/inf <!Tr 
k J R3 ) a 



T„l/2(A)—r—r+l 



\X\ 



+ - 



KV 



|V®A| 



We apply Theorem 12. 71 to the potential V(x) = — 1 and with h = v x l 2 . Notice that V[x) 
satisfies the conditions (12.161) and (12.171) with /i — 1, M — 1, r min — oo, z± — 1 and r x = 0. 
We get 



inf (Tr kfA)-^ 
ALL \x\ 

which proves (12.141) . 



~l |VM 2 }-7A 



P ~ TT + ?y 

9 



dpdq = 2S{k)+0{u £/2 ) 



□ 
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